function OpenCRTBP_MD_BaseFigure
%% Constants
G           = 6.674e-11 * ((1/1000)^3);
mass_Mars   = 6.39e23;
mass_Deimos = 1.8e15;
mass_Phobos = 10.8e15;
GM_Mars     = G*mass_Mars;
GM_Deimos   = G*mass_Deimos;
GM_Phobos   = G*mass_Phobos;
GM_sun      = 1.327124400*10^(11);     
%[km^3/s^2] Gravitational Parameters

% u_sm = GM_Mars / (GM_Mars+GM_sun);
u_md = GM_Deimos / (GM_Deimos+GM_Mars);
% u_mp = GM_Phobos / (GM_Phobos+GM_Mars);
%[] Mu constant
%    sm --> Sun-Mars system
%    em --> Earth-Moon system%%

% DU_sm = 235.34e6; 
% DU_mp = (9236.53+9517.58)/2;
DU_md = (23455.5+23470.9)/2;
%[km] Distance from Sun to Mars

% TU_sm = sqrt(DU_sm^3/GM_sun);
% TU_mp = sqrt(DU_mp^3/GM_Mars);
% TU_md = sqrt(DU_md^3/GM_Mars);
%[s] Time unit 
%[] Mu constant
%    se --> Sun-Earth system
%    em --> Earth-Moon system
rd = 6.2/DU_md;
rm = 3389/DU_md;
figure
hold on
[L1,L2,L3,L4,L5] = librationPoints(u_md);
Lse = [L1,L2,L3,L4,L5];
for ii = 1:5
    plot3(Lse(1,ii),Lse(2,ii),Lse(3,ii),'bo','MarkerFaceColor','b','MarkerSize',2)
end
[xs,ys,zs] = ellipsoid(-u_md,0,0,rm,rm,rm,50);
fv = stlread('Deimos.stl');
fv.vertices = (fv.vertices)/DU_md;
fv.vertices(:,1) = fv.vertices(:,1)+(1-u_md);
patch(fv,'FaceAlpha',1,'EdgeColor','none');
surface(xs,ys,zs,'FaceColor',[255 165 0]/255,'EdgeColor','None');
xlim([-1.3,1.3]);
ylim([-1.3,1.3]);
axis equal
set(gca,'FontSize',20) 

end